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A theoretical investigation of the trion formation process from free carriers in a single GaAs/ 
Ali-^Ga^As quantum well is presented. The mechanism for the formation process is provided by the 
interaction of the electrons and holes with phonons. The contributions from both the acoustic and 
optical phonons are considered. The dependence of both bi-molecular and tri-molecular formation 
rates on temperature is calculated. We demonstrate that they are equivalent for negatively and 
positively charged excitons. 
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In semiconductor wells, the photoluminescence spec- 
trum following the generation of electron hole plasma 
is dominated by an exciton lme-££££&. The for- 
mation of excitons in quantum wells has been ex- 
tensively investigated both experimentally-^i&ii&ji and 
theoreticall y 10 ' 11 ' 12 . Recently, it was experimentally 
showed that the exciton formation is strongly density 
and temperature dependent; it is a bimolecular process in 
which an electron and a hole are bound by Coulomb in- 
teraction with the emission of the appropriate phonon 8 . 
This experimental result unambiguously confirmed the 
theoretical prediction for luminescence spectrum is mod- 
ified: a charged exciton resonance appears below that of 
the exciton. 

In a previous publication^, we showed experimen- 
tally that the dynamics of exciton, trion and electron-hole 
plasma can be ruled by a simple rate equation model, in 
which we account for bimolecular formation of excitons 
from an electron-hole plasma, bimolecular formation of 
trions from excitons and free carriers and trimolecular 
formation from free carriers. Using only two fit param- 
eters, we were able to extract experimentally the depen- 
dence of both bi-molccular and tri-molecular formation 
coefficients on temperature. 

In this paper, we propose a theoretical derivation 
of these coefficients. We consider formation channels 
through which the formation channel and show that they 
correspond to experimental calculations. 

In Sec. I, we first reproduce the results of Piermaroc- 
chi et al.— on exciton formation. We then extend the 
formalism to the case of charged excitons; in Sec. II, 
we calculate the bi-molecular formation and in Sec. Ill 
the tri-molecular formation. In Sec. IV, we show how 
these formation rates can be used to calculate formation 
dynamics at thermodynamical equilibrium. 



The dynamics of the exciton formation is considered 
in the framework of the Boltzmann equation for a sys- 
tem containing free electrons, free holes, and excitons. 
The residual Coulomb interaction between the free car- 
riers is neglected, which is justified in the range of tem- 
peratures and densities considered—. In this work we 
focus on the exciton formation mechanism and do not 
discuss the relaxation of the three species within their 
respective bands, the electron-hole scattering, and radia- 
tive recombination. We denote the occupation numbers 
for electrons, holes, and excitons by f e (k e ), fh(kh), and 
fx(kx), respectively, where k e , kh, and kx are the in- 
plane momenta for electrons, holes, and excitons. For the 
bimolecular formation, following Piermarocchi et al.— , 
the scattering terms in the Boltzmann equation process 
reads 



d/ e (fc e ) 

At 



form 



E 

k x ,k h 



w ke ,k h ^k x fh(kh)f e (k e ), 



(1) 

where Wk B .k h ^k x represents the probability per unit 
time for a free electron and a free hole to bind to- 
gether and form an exciton. Free carriers thermalize very 
quickly in comparison to the exciton formation time, no- 
tably through fast carrier-carrier scattering 1 -.. It is thus 
assumed that during the evolution of the system, the free 
electrons and holes are thermalized at the same temper- 
ature T c . In the scattering term of Eq. |T]), we use for 
/ e (fc e ), and fh{kh) equilibrium Boltzmann distribution 
function at T c . Consequently, by summing Eq. §Q over 
fc e , we obtain an adiabatic equation for the evolution of 
the electronic density n e = j;J2k e fe(k e ) 



'Y^F(kx)n e n h = -Cn e n h 

kx 



(2) 



The coefficient C is the bimolecular formation coefficient, 
which depends on both T c and the lattice temperature T; 
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through the term 



F{kx) = 



2nh 2 
kuT r 



1 



k s ,k h 



m e m h S 

-(E e (h e )+E h (k h ))/k B T c /g\ 



where S denotes the QW surface area. 

Free carriers are coupled to the exciton by a contin- 
uum of phonon states (q,q z ) through a carrier-phonon 
interaction Hamiltonian T~l e /h-ph- A phonon can be 
emitted (+) or absorbed (-) in the formation process of 
the exciton. We calculate both case separatly using the 
Fermi's golden rule 



2tt 



w 



k,,k h 



E 



\(k X \ ® (n q ,q z ± l\H e/h - ph |«q,g 2 ) <B> \k e ) <8> \k h )\ 

x 6[E e {k e ) + E h (k h ) - E x (kx) T tiu ph (q,q z )], (4) 

with E e (k e ), Eh(kh) and Ex(kx) the energy disper- 
tion of the electrons, holes and excitons respectively, and 
ftw p h(q, q z ) the energy of the emitted (absorbed) phonon. 
We first build the bound and unbound electron-hole pair 
states \kx) and \k e ) ® \kh). 



A. Bound and unbound exciton states 

Let (r e ||,z e ) and (r h \\,z h ) be the electron and hole 
position vectors respectively and & x [(r e y , z e ). (r h \\ , z/J] 
the exciton wavefunction, where we have separated the 
coordinates in the QW plane (x-y) from the perpendicu- 
lar coordinates (z). Denoting the electron (hole) in-plane 
momenta k e (k^), we write the in-plane Fourier trans- 
form of this function 

^h e ,k h ( z ^> z h) =g I dr e\\dr hll <f> x [(r eh z e ),(r hh z h )} 



x e 



i(ka-r e »+k h -r h «) 



(5) 



where S denotes the QW surface area. Transforming to 
center-of-mass (CM) and relative coordinates in the QW 
plane — R\\ = a x r e \\ + (3xr h \\, r\\ = r e \\ - r h \\, where 
a x = m e /M x , fix = m h /M x and m e , m hl M x are 
the electron, hole and exciton in-plane effective mass — 
we can apply Bloch's theorem and decompose the exci- 
ton wavefunction into a free motion part e lkx related 
to the exciton in-plane momentum kx and an envelope 
function. To facilitate the calculation, we use an enve- 
lope function separable in z and rii , although it is strictly 
justifiable only for narrow well structures, 



(j) {r\\,z e ,z h ) = Xe{z e )xh{z h )(p x {r\\). 



(6) 



The confinement functions Xe(z e ) {Xh( z h}) is taken to be 
the wavefunction of an electron (hole) in the ground state 



of a finite square quantum welli^ 

A a cos (kl"^ z a ] 

Xa(z a ) = * 



for \z a \ < 



B a exp (-k h z (a) {\z a \ - L z /2)) for \z a \ > !f 

(7) 

with a — e,h and we use the simplest electron orbital 
function 



nX x 



whose in-plane Fourier transform is given by 



'87rAx z 



[1 + (X x k) 2 ] 



21-3/2 



(8) 



(9) 



The variational parameter Xx is associated with the Bohr 
radius of the exciton in the QW. Eq. © can be rewritten 
as 



1 



$k e ,k h ( z e>Zh) =- / di?||dr||0 (r\\,z e ,z h ) 



x e ik x -R\\ e -i[R\y(k e +k h )+rw(l3xke-a x kh,)} 



(10) 



where 4> x (z e ,z h ) = Xe{z e )Xh(zh)<P k x is tne in-plane 
Fourier transform of the exciton envelope function. We 
can now construct the state of a single exciton with an 
in-plane momentum kx in the Fermionic Hilbert space of 
electron- hole pairs. It is the superposition of wavefunc- 
tions (dU]) with all electron momenta k e and all electron 
z e and hole Zh coordinates, given by 



\kx) 



k c 

x c T 



dz e dz h (t) x * kx+kc (z e ,z h ) 



|0>, 



(11) 



where c£ z (d kx z ) is the electron (hole) creation oper- 
ator with in-plane momentum kx and z e (z^) coordinate. 

Similarly, we choose plane waves for the free carriers. 
Thus the in-plane Fourier transform ip k {z) of the carrier 
wavefunction takes the simple form 

0?W=XaW, a = e,h. (12) 
The unbound electron-hole pair then reads 

\ke,k h ) = J dz e dz h ^(z e )^l(z h )c{ etZ S khiZh |0>. 

(13) 

B. Carrier-phonon interaction Hamiltonian. 

We write the interaction Hamiltonian for a coupled 
electron-phonon system in the notation of the second 
quantization 16 

U a - P h = ^2 V q,q* Qa(q,qz) (14) 
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where a q ^ is the phonon creation operator. The 
electron density operator Q a (r\\,z) and its counterpart 
in Fourier space Q a {q, q z ) are expressed on the basis 
{<f > %(r ll ,z ) s) = e ih ^5(z)Us)} 



Qh( r b z ) 



k,k' 

Ql{q,q z ) = Y J [ dz e-^%l q 

k 



k.z 

i{k — fc')-r|| 



k.k' 



u k+q,z u k,z- 



(15) 

(16) 
(17) 

(18) 



where eo is the static dielectric constant and e x is the 
high frequency dielectric constant. We use the notation 
a a for the deformation-potential constant (assumed to be 
associated with a non-degenerate conduction or valence 
band), po for the density of the crystal, e for the charge 
of the electron and V for the volume of the sample. We 
follow Einstein interpolation scheme, so that the disper- 



sion is merely w 



LA 



,LO _ 



= ujlo for LO phonons, 



ql for LA phonons and 
standing for the Debye 



sound velocity and ujlo for the reststrahl frequency. 



Spin states Cx( s ) = ( cr l s ) have been intruduced for their 
will be necessary when we treat the trion formation. 

Only longitudinal acoustical (LA) and longitudinal 
optical phonons (LO) couple significantly to careers. We 
express the coupling vertex functions V qq ^ for both cou- 
pling 



ya (LA) ■ 



/ ft(!g| 2 + gf) ; 

2p Vuj q ^ z 



(19) 



ya(LO) = 

V (\q\ 2 + q 2 z )V \e e 



Matrix element calculation. 



The matrix elements in Eq. (T4)) are calculated, mak- 
ing use of Eq. (JTIJ, (TT3]) and (33]) 



{k x ,n q ^ qz ± 1| 7Y c _ p h + Hh-ph \k e , ku, n q ^) 

= Y2J2 f dz e dz h dz edz' h dzipX kx+K (n q . qz ± 1 1 a ±qAz al -_- J n qAz ) X * e (z e )Xe {K)Xh ( z h )Xh (4 )e~ lqzZ 

q,q z 

x{vi 4z (0\c. K , z ^ k+q ^c^ 
Applying operators on the ground state 

(0\C-k' e , Z ' e C k ^ Ze \0) = 5k c .-k' a S(z e - z' e ), 

zh), (0|4x+^,<4 h , z J°> = Sk h ,k x +k' c 8(z h - z' h ), 



we obtain 

(k x ,n q ^ z ± 1| Hc-ph + Tih-ph \k e ,kh, n q . qz ) 

= \j n i,q, +1 ± 1 5 ±qM e +k h -k x {vi qz tp X p xkx+kh I e {q z ) - V^ q ^ xhx _ he I h {q z )^, (21) 

where the integrals, in the orthogonal direction are given by 

I a (q z )= [ dz a \ Xa (z a )\ 2 e lq ^, a = {e,h}. (22) 



(0\c-k' e , Z ' s c k+q z Ck, z c k<! Z<! |0) = 8-k^.k+q 5{z - z' e ) 5k.k c S(z- z e ), 
(0\dk x +k' c , z 'Ji + q,Jk, z d{ htZh \0) = 5k x +k' e ,k+ q 5(z - z' h ) 8 k ,k h 8(z - 
(n q<qz ± l\a q , qz + _-Jn g , 9 J = Jn q . qz + \±\ &F9,9&F««,g*> 
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Finally, if we choose the bound and unbound electron-hole pairs dispersion to be parabolic, the probability transi- 
tion (0J reads 



2tt 



w 



k,,kj.— >fej 



X d 



2m e 2ra} l 



xl-x -kjh{<lz)\ 
hi. 



2M 



5±q,ke+k h -k x - (23) 



D. LA phonons assisted formation 



and reexpressing the Dirac distribution as 



Considering that the sample volume V — L Z S 
is macroscopic the sum over the orthogonal phonon 
wavevectors may be replaced by the integral 



E 



2tt 



(24) 



2 1,2 



2m e 



2rrih 



\?iv s q z 



-p {s\q z -q z Hke,k h ,kx))} + 5 \q z + qf\k e , k h , k X j\ } 

2 L L J L J J 



± 



h 2 k 2 

2m„ 



2m h 



h 2 k\ 
2M 



(25) 



with 



q z 0) (k e , k h ,k x ) 



1 



h 2 k 2 h 2 k 2 

h 2 v 2 \ 2m e 2mh 
I 



h 2 k\ 



q 2 , 



(26) 



makes the integration (j2~4"|) trivial for LA phonons: 



E. LO phonons assisted formation 



4ttL ? 



Wk a .k h ^k 



q'+q^ 



h 2ir2pVv s \hv s qi 0) \ 
aef0 x k x -k,Je(q z O) ) - a h (Pa x k x -kJh(q z 0) ) 

h 2 k 2 



In (o) 



2m„ 



2m h 



+ n (o) 



2m P 



h 2 k 2 
2m h 



h 2 k 2 

E b £ 

2M 



h 2 k 2 x 

2M 



(27) 



This expression already includes the sum over absorbed 
and emitted phonon contributions. The phonon in-plane 
momentum needs to be substituted by q = k e + k b — kx- 



In the case of interaction with LO phonons, Eq. 
becomes 



Wk a ,k h ^k x 



2ir hoj LO e 2 {l/e 00 - l/e ) 



V(q 2 + q 2 ) 



X IVfSxkx-kJeiqz) - ¥ ax k x -kjh{qz)\ 

'h 2 k 2 



x S 



2m P 



2m h 



E b 



h 2 k 2 x 

2M 



huj LO 



(28) 



where we dropped the phonon absorption part, which is 
negligible up to room temperature. For the calculation 
of the exciton formation coefficient C, it is convenient to 
rewrite the Dirac distribution as 



h 2 k 2 



h 2 k 2 



2m e 



H 2 k 2 

E b ^ 

2M 



2nih 

^{S[k x -kf(k e ,k h )) 



LO 



k x + k 



(k e ,k h ) | 



(29) 
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FIG. 1: The exciton formation coefficient C as a function 
of the carrier temperature T c , at a fixed lattice temperature 
Tj = 10 K. Other parameters are given in the text. 



Parameter 


Symbol Value 


Unit 


Band gap energy 


E g 


1519 


meV 


Electron effective mass 


T7le 


0.08 


mo 


Heavy hole effective mass 


m hh 


0.17 


m 


LO phonon energy 


hulLO 


36 


meV 


Static dielectric constant 


CO 


12.85 




High frequency dielectric constant 


too 


10.88 




Crystal density 


P 


5.3162 


g-cm -3 


Sound velocity 


v s 


4726.5 


m-s -1 


Conduction band deformation potential 


a e 


-7.0 


eV 


Valence band deformation potential 


ah 


3.5 


eV 


Exciton binding energy 


E T 


6.5 


meV 


Trion binding energy 


Et 


1.77 


meV 


Exciton Bohr radius 


\x 


11.0 


nm 


Trion variational parameter #1 {X~) 


At 


15 


nm 


Trion variational parameter #2 (X~) 


Ay 


30.0 


nm 


Trion variational parameter #1 (X + ) 


At 


16 


nm 


Trion variational parameter #2 (X + ) 


Ay 


25.0 


nm 



TABLE I: GaAs Material Parameters 



with 



>2Mx (h 2 kl 



h 2 



2m P 



2m h 



Eh — huj 



Numerical Results. 



LO 



(30) 



where Fk ht k e ^k T represents the probability per unit time 
for a bound electron-hole pair and a free electron to bind 
together and form an exciton. We assume that bound 
and unbound carriers are thermalized and use Boltzmann 
distribution function fx{kx) and f e (k e ) for exciton and 
electron population, respectively. By summing Eq. (|31[) 
over fe e , we obtain an adiabatic equation for the evolution 
of the electron density 



If we change all the sum in Eq. H El [37] and [25] into 
integrals, the bimolecular formation coefficient C can be 
numerically calculated by Monte Carlo integration. In 
Fig. [U we report C as a function of 1 /T c for a fixed lattice 
temperature 2] = 5 K, for a GaAs QW of 80 A. The two 
contributions from the acoustic and optical phonons are 
shown separately. The acoustical phonon dominates for 
temperatures smaller than 40 K and does not depend 
on T;. We see that these results perfectly match those 
published by Piermarocchi et al.— . We show in Table U 
the numerical value of the different parameters entering 
in the calculation. 



dn 



n x n e 



-A 2 n x n e 



(32) 



The coefficient is the bimolecular formation coeffi- 
cient, which depends on both T c and the lattice temper- 
ature Ti through the term 



F(k T ) 



2i:h 2 



k B T c J m x m e S 
x e 



(34) 



k x ,k e 

(E x (k x )+E.(k B ))/k B T, 



III. BIMOLECULAR FORMATION OF TRIONS 

We now extend our formalism to the derivation of 
the bimolecular formation of trions. We restrict ourselves 
to negatively charged excitons and will give at the end of 
this work some indication on how to retrieve their posi- 
tive counterpart. 

We write the scattering term in the Boltzmann equa- 
tion process for the bimolecular formation of trions: 



d/ e (fc e ) 
At 



form 



= - F k X ,k s ^k T fx(kx)fe(k e ), 



(31) 



We calculate the formation rate using Fermi's golden rule 



2tt 



w 



k x -k e — >kr 



E 



\(k T \ ® (n q ,^ ± 1| U e/h ^ ph |n g , 9a ) ® \k x ) <8> |fe e )| 
x S[E x (k x ) + E e {k e ) - E T {k T ) T ?Wph(q, Qzj] , (35) 

with E x (k x ), E e (k e ) and Ex^t) the energy disper- 
tion of the electrons, holes and excitons respectively and 
ftw p h(q, q z ) the energy of the emitted (absorbed) phonon. 
We first build the bound and unbound electron-hole pair 
states \kr) and \k x ) (g) |fc e ). 
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A. Trion state. 

The two electrons and the hole are positioned at 
{ri\\,Zi), (*"2||,Z2) and { r h\\, z h) respectively, while the 
center-of-mass (CM) and relative coordinates in the QW 
plane are now given by U|| = axi^iw + f 2 ||) + PrfhW and 
nh|| = rq - r h\\ {i = 1,2); a T = m e /M T , (3 T = m h /M T 
and Mt is the trion mass. We consider a simple two 
parameter Chandrasekhar-type trial envelope function 
that was successively used to calculate trion-electron 
scattering^: 

T (»"lh|b r 2h\\, z l,Z 2 ,Z h ) 
= NTXe {Z\ )Xe {z2)Xh {z h ) 

[<p XT (r lh )<p x 'T(r 2h ) ±<p X ' T (r lh )v XT (r 2h )], (36) 

where the + (— ) sign applies to the singlet (triplet) spin 
configuration and the trion wavefunction normalization 
factor is given by 



Mr 



1 



with 



V2(l±« 2 ) 



4AA' 



(A + A') 2 ' 

Its in-plane Fourier transform reads 



(37) 



(38) 



ki,k 2 ,k h ( Z ll Z 2) Z h) 



Sk T -ki-k 2 -k h t f > aTk T -k 1 ,a T k T -k 2 ( z 'L> ^2, Zh) (39) 



and the state of a single trion with an in-plane CM mo- 
mentum kx is constructed similarly to that of a single 
exciton 

| fc f ) = j dZ l dZ 2d Z I.CfcT + ltl,«TfcT+fc2( Z l' Z 2' 

SI ,S2 

x s 2 ) c s X Zl c s X,A T +k 1+k2 ,, h |0> , (40) 



Zh ) 



where we have added the spin index to the electron cre- 
ation operator and introduced £s(si, s 2 ) = (<S|si, s 2 ) the 
projection of a generic spin configuration of two electrons 
on the singlet spin configuration. The in-plane Fourier 
transform of the trion singlet wavefunction in Eq. (|39p 
is given by 



4>ki,k 2 (Zl ,Z 2 ,Z h ) =Af T Xe(zi)Xe {z 2 )Xh (zh) 



(41) 



where <p k ^ has already been defined in Eq. [9l The 
exciton-free electron state is given by 

\k X ,k e >)=Y^ J dz e dz h dz' e <t>a x kx+ke( Z e>Zh)'>l>h e ,( Z e') 



x cL fca ^4x+*.,. h 4^,,.J0>. (42) 
B. Matrix element calculation. 



We calculate the matrix elements in (1331) 



(fcf;s/J ® {n q ^ ± l\H e / h -ph \ n q,g z +1) ® |fex;si;s ft ) ® |fc 2 ;s 2 ) 



X! X! X! X! J dz dzidz 2 dz hdz^dz^dz'h \/n q . qz + 1 



(°l cic 2 4c 4 44 10) (0| d r 4 |0) - , 



where we simplified the operators index using the following scheme 

1 = ( — z i, s^) 2 = (~k 2 , z 2 , s 2 ) 

4 = (k,z,s) 5 = (-fei,zi,si) 

7 = (&t + fe'i + fc 2 , z ' h , s' h ) 8 = (fcx + ki,Zh, Sh) 



(0|d 7 444l°> • (43) 



3 = (fc + q, z, s) 
6 = (k 2 ,z 2 ,s 2 ) 



Using the electron, holes anti-commutation relations, the Fermi vacuum expectation value of the operators read 

(0| C^C^cJcJ 1°) = ^13(<5 2 5^46 - S 26 S 45 ) - <5 2 3(<5l5^46 ~ SlB^As) (44) 



<0| d 7 d\ |0) = s 78 
)| c^clcl |0) = S 25 6 16 - S 15 6 26 
I d 7 d\d 4 d\ |0) = 63764$ 



(45) 
(46) 
(47) 
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Using the later results in Eq. (|43|) gives 

(fef ; sJJ O (n g ,g a ± 1| H e / h -p h \n g ,q z ) <8 |fex; si; s;i) ® Ifo; s 2 ) 



= si) - £(«i,* a )] S Sh ,s' h ^n q<qz + 1/2 ± 1/2 £ ^ xfex+fel 

fei 

X {y r q,q^^ e (l z )l- ( t'a T kT-k2,aTkT+ki+q + ^a T fc T -fc 2 +g,aTl!T+feJ _ ^q,q z ^h{Q_z)<t>a T k T ~k- 2 ,a T k T +k l ^^qMx+k 2 -k 2 

Finally, for parabolic electron, exciton and trion dispersion, the probability transition (|35p is 



fex ,k 2 —>k T 



q,q* 



^ ] ^,g2^ e fe)[0a T fcT-fc2,ctTfcT+fcl+q 0a:rfeT— fc2+g,arfer+fei] ^ (.Q z ) C t ) aTk T -k 2 ,a T k T +ki 

fel 

,21.2 



X (5 



2M 



x 



2m P 



2M 7 



=F hw ph (q,q z ) 



5±q. 



q,kx+k 2 —kj 



(48) 



(49) 



where we averaged over the initial electron spin states and e xciton angular momentum states. 



IV. TRIMOLECULAR FORMATION OF 
TRIONS 

Again, we write the scattering term in the Boltz- 
mann equation process 



fc e , we obtain an adiabatic equation for the evolution of 
the electron density 



~d7 



^2F(k T )n 2 e n h = -A 3 n 2 e n h . (51) 



d/e(fc e ) 
dt 



The coefficient A 3 is the trimolecular formation coeffi- 



= — Ffc e .fc^fc h ^fc T / e (fc e )/ e (fcg)/^(fc/j),cient, which depends on both T c and the lattice temper- 

form hr,K 



(50) 

where Fk e ,k' ,k h ^>k T represents the probability per unit 
time for two free electrons and one free hole to bind to- 
gether and form an exciton. We assume that bound and 
unbound carriers are thermalized and use Boltzmann dis- 
tribution function fx(k e ) and fh(kh) fo r exciton and elec- 



ature TJ through the term 

-2x3 



F{k T ) 



(2nt 



-(Ee(k e )+E e (k' e )+E h (k h ))/kBT a 



w k c ,k e ,k'^ki 



X e 



(52) 



tron population, respectively. By summing Eq. (|50[) over We calculate the formation rate using Fermi's golden rule 



2tt 

K,k h ^k T = jE K fcT ' ® ( n 1^ ± MKe/h-ph \nq,qz) ® |*!e) ® \K) ® 



5[£ e (fce) + E e {k' e ) + E h (k h ) - E T {k T ) T fiw ph («i, a,)] , (53) 



with Ex(kx), E e (k e ) and Erikr) the energy dispertion of the electrons, holes and excitons respectively, and 
hw p h{q, q z ) the energy of the emitted (absorbed) phonon. Finally, for parabolic electron, exciton and trion dispersion, 
the probability transition [53] becomes 



2tt 



,ki,kh—>k'i 



~^^( n q-qz + 2^2' 



q,q* 



Vq,qzIz( < lz)[4'a T kT-k2,a T kT+k 1 +q + ^arkr -k 2 +q,u T k T +k 1 ] ^q,9z ^ fes ) < t > u T k T ~k 2 ,a T k T +ki 



x (5 



rft 2 fc e 2 
2m P 



2^2 



h z k 



2m e 



2m h 



h 2 k 2 



S±q,k x +k 2 -k T , (54) 
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where we averaged over the initial electron spin states and e xciton angular momentum states. 



A. Numerical results 
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FIG. 2: The bi- and tri- molecular trion formation coefficients 
A% and ^3 as a function of the inverse carrier temperature as 
calculated from our model for LA phonon assisted formation. 
The red curve indicates X + formation and the blue curve X~ 
formation. 

In Fig. [2j we represent the results of our numeri- 
cal calculation for LA phonon assisted formation. We 
stress the fact that the results for LO phonons are or- 
ders of magnitude smaller and consequently negligible 
for bi- and tri- molecular processes. This shows that the 
exclusive formation mechanism for trions is governed by 
LA phonon interaction. We predict a decrease of the bi- 
or tri-molecular formation for raising temperatures. We 
also demonstrate that the bi and tri-molecular formation 
coefficients for negatively and positively charged excitons 
are approximativcly equal (A 2 = A 2 ) and (A3" = A3" ) . 

A more accurate calculation should rely on a bet- 
ter trion wavefunction. The Chandrashekar variational 
function is most likely to simple to yield quantitative re- 
sults. We note however that electron-trion scattering 18 
will contribute to ionize trions. This effect will dramat- 
ically increase at high temperatures. We consequently 
propose that the electron-trion scattering reduce the ex- 
perimental values for the trion formation coefficients. 



V. FORMATION RATES CLOSE TO 
EQUILIBRIUM 

In this Section, we propose to derive formation 
rate for electrons, holes, excitons and trions assuming 
thermodynamical equilibrium. In this set of equations, 
we neglect both biexciton channels and Auger channels 
because experiments on undoped samples demonstrate 
that, at the densities considered in the present work, 
these channels are not significant 8 -. 

The dynamics of a plasma of electrons containing 
electrons (e), holes (h), excitons (X) and trions (X + and 
A~) is governed by the following five channels: 



e + h <-» X 

X + e^X- X + h^X+ 

2e + h^X~ 2h + e^X + 

The different formation rates for these populations 



read 



F x = Cnp - jCK x X, 



= A 2 -Xn-A 2 -if 2 -X- 
= Ag~nnp — Ag" A^X - , 
= A+Xp-A+A+X+, 
= A+npp-A+A+X+, 



(55) 
(56) 
(57) 
(58) 
(59) 



where C, A 2 and A3 are respectively the exciton, 
trion bimolecular and trion trimolecular formation rate 
calculated in this article; E^x is the exciton binding 
energy; Ax, K 2 , Kg the equilibrium coefficients. For 
a 2D system, they can be derived from the Boltzmann 
distribution 



•) . / „-hbar 2 k 2 /2m a k B T _ 

Z9a * 2n 



m a k B T , KbT 
9a 27rfi 2 6 



(60) 

where a = e, h, X, T, the factor g a = 2 is the spin degen- 
eracy of the electron, hole, exction and trion in the non 
degenerate regime, where Boltzmann statistics applies. 

Using that fact that the chemical potential of the ex- 
citon nx is related to the chemical potential of electrons 
He and holes Hh as fix — Me+M/i+^x (Ex is here defined 
as minus the exciton binding energy), one immediately 
obtains the Saha relation for the exciton density 

n? = Kx[T) = g^rn^k^T /kBT 
X gx mx 2-nn 

Similarly, the chemical potential of the trions is 
fix = 2/i e + fj,h + Et (Et is here defined as minus the 
trion binding energy), so that Saha equations for trion 
bi-molecular formation are 

Xn - Z—(T-\ - ffeffx rn x m e k B T „- { E T -E x )/k B T 

(62) 

Xp - K+(T\ - 9 fa* m xm h k B T ^_ (ET _ Ex)/kBT 

(63) 



X+ 



Finally, the following set of equation is infered for 
tri-molecular formation: 



n 2 h 

IF 



gr mx- \2irn z ) 



— = K + (T) = - 9eg ft me7n h ( k ' BT 
X+ 3 [ ' g T m x + \2irh 2 



-E T /k B T 



(65) 
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Such a set of equations, together with the equilib- 
rium densities, allow to compute the dynamics of the dif- 
ferent populations after non-resonant optical excitation. 
This may apply both to undoped^ as well as to doped 
quantum wells. The results in the case of a sample doped 
with electrons will be detailed in another publication 13 . 
The experiments show that indeed, when the density of 
electrons is sufficient and at excitation densities of the 
order of 10 10 cm -2 and above, the trimolecular formation 
process of trions has to be taken into account to properly 
reproduce the observed dynamics. 

VI. CONCLUSION 

In this paper, we have derived the appropriate model 
for computing the rates for exciton and trion formation. 
For the case of trions, we have derived the equations for 
both the bi- and tri-molecular phonon-assisted formation 
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